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Abstract 

The effect of temperature inhomogencity on the periods, their ratios (funda- 
mental vs. first overtone), and the damping times of the standing slow modes in 
gravitationally stratified solar coronal loops are studied. The effects of optically 
thin radiation, compressive viscosity, and thermal conduction are considered. 
The linearized one-dimensional magnetohydrodynamic (MHD) equations (under 
low-/3 condition) were reduced to a fourth-order ordinary differential equation 
for the perturbed velocity. The numerical results indicate that the periods of 
non-isothermal loops (i.e. temperature increases from the loop base to apex) 
are smaller compared to those of isothermal loops. In the presence of radiation, 
viscosity, and thermal conduction, an increase in the temperature gradient is 
followed by a monotonic decrease in the periods (compared with the isothermal 
case), while the period ratio turns out to be a sensitive function of the gradient 
of the temperature and the loop lengths. We verify that radiative dissipation 
is not a main cooling mechanism of both isothermal and non-isothermal hot 
coronal loops and has a small effect on the periods. Thermal conduction and 
compressive viscosity are primary mechanisms in the damping of slow modes 
of the hot coronal loops. The periods and damping times in the presence of 
compressive viscosity and/or thermal conduction dissipation are consistent with 
the observed data in specific cases. By tuning the dissipation parameters, the 
periods and the damping times could be made consistent with the observations 
in more general cases. 

Keywords: Sun, corona; Sun, magnetic fields; Sun, oscillations 
1. Introduction 

Through SOHO and TRACE observations, propagating slow magnetoacoustic 
waves have been found in coronal plumes, loop footpoints, above sunspots and 
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even non-sunspot regions (e.g., Ofman et al, 1997; DeForest and Gurman, 
1998; Berghmans and Clette, 1999; De Moortel et al, 2000; OShea et al, 2002; 
Brynildsen et al, 2002; Marsh et al., 2003). Standing longitudinal slow waves 
with strong damping and large Doppler-shift oscillations have been detected 
in hot postflare loops recorded by SOHO/SUMER. In the cooler loops, these 
waves were observed by the EUV imaging spectrometer Hinode/FAS (Srivastava 
and Dwived, 2010). These oscillations have a phase shift of about one-quarter 
period between velocity and intensity. The periods and damping times of the 
standing slow waves are in the ranges of 8.6-32.3 min and 3.1-42.3 min, re- 
spectively (Kliem et al, 2002; Wang et al, 2002a, 2002b, 2003, 2005; Banerjee 
et al., 2007; Erdelyi et al., 2008). The effect of energy dissipation on the slow 
waves through thermal conduction, compressive viscosity, radiative cooling, and 
heating on the periods, period ratio, and damping times has formed the focus 
of recent studies. Ofman and Wang (2002), for instance, studied the oscillations 
and damping of standing slow modes in the isothermal loops. They found that 
thermal conduction is the dominant dissipation mechanism. Taking into account 
the effects of thermal conduction and compressive viscosity, De Moortel and 
Hood (2003) investigated both propagating and standing slow magnetoacoustic 
waves in homogeneous corona loops. They found a minimum damping time that 
can be obtained by thermal conduction alone. However, for stronger dissipation 
an additional mechanism such as viscosity has to be added. Sigalotti et al. 
(2007) studied the dissipation of standing slow modes in hot, isothermal loops by 
integrating the effects of gravitational stratification, thermal conduction, com- 
pressive viscosity, radiative cooling, and heating in their study. They concluded 
that thermal conduction and compressive viscosity are the main sources of the 
wave damping. Pandey and Dwivedi (2006) have shown that separate effects of 
thermal conduction and viscosity in isothermal loops are not sufficient to explain 
the observed damping times of the oscillations. Only with the combined effect 
of thermal conduction and viscosity the results are consistent with the obser- 
vations. Moreover, Taroyan et al. (2005) investigated the effect of temperature 
inhomogencity on the dissipation of standing slow waves through considering 
thermal conduction and optically thin radiative losses. They found that the 
damping time of the isothermal loops is proportional to the wave period and the 
oscillations are rapidly damped mainly by thermal conduction. 

The ratio between the period of the fundamental mode (p\) and the first 
overtone (P2), Pi/(2p2)j for both the fast kink mode and the slow longitudinal 
mode is a useful scismological indicator. Departure of this ratio from unity is a 
consequence of the density, pressure, temperature, magnetic field (both radial 
and longitudinal structures), heating functions and the dissipation mechanisms 
( Diaz et al., 2006; Dymova and Rudcrman, 2006; Donnelly et al., 2006; Erdelyi 
and Vcrth, 2007; Safari et al, 2007; Vcrth et al, 2007; Rudcrman et al, 2008; 
McEwan et al, 2008; Andries et al, 2009; Fathalian and Safari, 2010). The 
results of the numerical modeling of the oscillations allow a comparison between 
the model and the observational results (McEwan et al, 2006; Abedini and 
Safari, 2011). Macnamara and Roberts (2010) studied the effects of thermal con- 
duction and compressive viscosity on the slow mode of the isothermal loops, and 
concluded that the effect of thermal conduction on the period ratio is negligible. 
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In line with the above-mentioned studies, in this study, we will investigate 
the dissipation of the standing slow MHD modes in the hot coronal loops taking 
into account gravitational stratification, inhomogcneity of temperature, thermal 
conduction, compressive viscosity, heating, and optically thin radiative losses. 

For such a purpose, the linearized one dimensional magnetohydrodynamic 
(MHD) equations are reduced to a fourth-order differential equation for the 
perturbed velocity. The paper is organized as follows. In Section 2, we present 
a brief description of the models and equations. In Section 3, the differential 
equation is solved numerically based on the boundary value problem solver (finite 
difference code) for different cases. Finally, conclusions are drawn in Section 4. 



2. Description of the Model and Equations 

Coronal loops are modeled as a magnetic flux tube of length 2L with a shape 
of half-circle and symmetrical around the apex (z = 0), with a constant mag- 
netic field along the loop length. The footpoints are fixed in the chromosphere. 
The equilibrium plasma flow is ignored, vo = 0. The projection of the gravity 
acceleration on the tangent to the loop (half-circle) is equal to gs'mirz/2L. The 
equilibrium temperature is considered as a function of z, 



T (z)=T (L) 



1 



TTZ 

(A-l)cos — 



- L < z < L, A = 



To(0) 



(1) 



where To(0) and Tq(L) are the temperatures at the apex and at the base, respec- 
tively. Under the ideal gas law as well as in the low-/? plasma, the equilibrium 
density and pressure are given by 



Po(z) = pa{L) |l + (A - 1) cos — 
Po (z)= Po (L) [l + (A-l)cos — 



, a T = 



tt(A-I) 



(2) 



where N = L/Aq(L) and the pressure scale height is defined by Ao = po/pog. 
The energy terms, the constant heating rate (H = Hq), the conductive heating 
(-Be), the radiative cooling (Er), and the compressive viscous heating (E v ) all 
are given by 



oz oz 



XP 2 T a , E v 



4 dv 2 
3 oz 



(3) 



in which k\\ = I0" 11 T 5 / 2 (Wm- 1 K- 1 ) and 77 = KT 17 T 5 / 2 (kg m^s -1 ) arc the 
coefficient of thermal conduction along the magnetic field and the coefficient of 
compressive viscosity, respectively (Braginskii, 1965; Hildner, 1974; Priest, 1982). 
The parameters a and x are a = —1 and \ = 10 24 Wm 3 kg _2 K, respectively, for 
the temperature range of 0.8 — 10 MK (Sigalotti et ai, 2007). In the equilibrium 



state, we get H + E ( 



oc 



Eh 



Orj 



-Eqr . For the convenience of our analysis, the 
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dimensionless parameters are defined as 

t ., - z - v - T p 

t = E/T s , Z = — , V = ; -, 1 = ; -, V = ; -, 

' ' V c s (L)' T (L)' 1 po(L)' 

v k - k \\ ^_ (i-i)t X pI{l)t-\l) 

1 _ m (L)r s (j - l)k oll (L)T (L)p (L) 

6 Re p Q (L)L^ 1 2 pI{L)t s 

where c 2 (L) = jPo(L)/ p (L) is the adiabatic sound speed, r s = L/c s is the 
sound crossing time from the loop footpoints to the apex, and Re is the Reynolds 
number. The dimensionless parameters d, r, and e are called the thermal ratio, 
radiation ratio, and the compressive viscosity ratio, respectively {e.g., De Moor- 
tel and Hood, 2003, 2004). In hot coronal plasma conditions, the values of r, e, 
d and t s are tabulated in Table 1. 



P 

Po(LY 

(4) 



Table 1. Dimensionless ratios and the sound crossing time as functions of 
the loop length (L) and the apex temperature (Tq(L)). 



T (L) [MK] 


L [Mm] 


r x 


IO 3 


e 


d 


r s [min] 


6 


20- 


-220 


0.05- 


-0.53 


0.02-0.21 


0.54-6.0 


0.88-9.64 


10 


20- 


-220 


0. Ol- 


-0.11 


0.08-0.83 


2.18-24.0 


0.69-7.65 



In the low-/? condition, the linearized non-ideal MHD equations are 
dpi . dvi dp 

-^ + Po a7 + aT 1 = ' (5) 

dv! dpx , iVsin(f ) Pl Ad 2 v x 
ot ryoz 7 3 oz z 

Pi _ pi T\ 
Po Po J-o 
dpi dp dvi 

-2rE^ + rE 0R ?± + d±(ko^ + k.A) . (7) 
po Jo oz oz oz J 

As we note here, all the quantities in Equations (5)-(7) are dimensionless (bars 
on dimensionless quantities have been dropped). The perturbed quantities are 
assumed to be of the form fi(z, t) = fi(z) exp(z4). After some algebra, Equations 
(5)-(7) reduce to a fourth-order differential equation as 

a, , ,d 4 vi ,d 3 vi „, , ,d 2 v\ 

A(z, r, e, d, v )~^ + B(z, r, e, d, v)-t-$ + C{z, r, e, d, v )~^ 
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+V(z,r,e,d,v)-^+£(z,r,e,d,v) = 0. (8) 

Coefficients A, C, T>, and £ are given in the Appendix. 
The boundary conditions are 

• The footpoints (z = ±1) are expected to be the nodes, i.e., f|footpoints = 0. 

• Near the apex, z — !> 0, all the coefficients in Equations (5)-(7) are constants 
and the solutions are expected to be v(z,t) ~ cxp[— i(wt — kz)\. 

In the remainder of this paper, Equation (8) is solved for five various cases by 
imposing the appropriate boundary conditions. 



3. Solutions to Equation (8) 

In all of the following cases, the loops arc considered to be gravitationally 
stratified. 

3.1. Case 1: d = r = e = 0, N ^0 

In the absence of all dissipation terms, Equation (8) reduces to a second or- 
der ordinary differential equation. The resultant equation is solved numerically. 
In Figure 1, the fundamental and the first overtone modes and their ratios, 
Pi/(2p2), are plotted versus L. As shown in the figure, in the case of adiabatic 
and isothermal loops (A — > 1), the fundamental and first overtone periods are 
in the ranges of 3-32 and 1-18 min, respectively. The results are well consistent 
with the previous studies (Curdt et al, 2002; McEwan et al., 2006; Sigalotti et 
al, 2007; Abedini and Safari, 2011). 

We see that, by increasing A from 1 to 2 the periods and their ratio decrease; 
yet, the changes in the fundamental periods are larger than those of the overtone 
periods. For typical loop lengths of 20-200 Mm, after adding the temperature 
gradient (A = 2) the changes in the periods is found to be about a few minutes 
(<10). 

3.2. Case 2: d = e = 0, r ^ 0, N ^ 

In the presence of optically thin radiation (i.e., e = d = 0,r =/= 0), similar to 
the previous case, we have a second order differential equation. The numerical 
solutions to the resultant equation are plotted in Figures 2 and 3. By comparing 
Figures 1 and 2, we find that the fundamental and the first overtone modes have 
changed about one percent. Therefore we conclude that the radiation does not 
change the periods and their ratios significantly. In Figure 3, the damping time 
and damping quality (damping per period, T<j/p = Re(w)/27rlm(w)) are plotted 
versus L for different Tb aS e and A. As shown in the figures, the damping time is 
of the order of 10 6 min and the damping quality is of the order of 10 5 . This case 
lies in the weak damping regime (r^/p > 2). By increasing A from 1 to 2, the 
damping time and damping quality increase. We arrive at the conclusion that, 
in the hot coronal loop, radiation is not a primary mechanism for the damping 
of the slow modes. 



S0LA1911_revl.tex; 6/03/2013; 17:42; p. 5 



Abedini et al. 



3.3. Case 3: r = d = 0, e 0, N ^ 

Next, the effects of compressive viscosity on the oscillations and damping of the 
loops are investigated. The numerical solutions are shown in Figures 4 and 5. 
By comparing Figures 1 and 4 we see the following: Compressive viscosity intro- 
duces a cutoff frequency on the oscillations as a function of both inhomogeneity 
parameter A and loop length L (see, e.g., Sigalotti et al., 2007). The effect of 
compressive viscosity on the periods for a typical loop length of 20-200 Mm is a 
few percent (< 2%) increase compared to case 1. By increasing inhomogeneity in 
the temperature (A) and decreasing the loop length, the ratio (pi/(2p2)) deviates 
from unity. At higher temperatures (> 6 MK) and non-isothermal loops, viscous 
dissipation is one of the primary damping mechanisms. Indeed, viscosity is one 
of the main factors for the strong damping in the short and non-isothermal loops 
and is consistent with the observational data. 

In the presence of the temperature gradient, the damping of the first overtone 
is stronger than that of the fundamental mode. Unlike radiation (case 2), if A 
is increased from 1 to 2 in the presence of viscosity, the damping time and 
damping quality decrease. We conclude that, for non-isothermal loops in the 
presence of viscosity and gravity, the computed damping times are consistent 
with the SUMER observations (Wang et al, 2003, 2005). 

3.4. Case 4: d ^ 0, r = e = 0, N ^ 

In the presence of thermal conduction, the periods decrease slightly if the tem- 
perature increases (Figure 6). However, by increasing inhomogeneity in the 
temperature (A from 1 to 2), the periods and their ratio decrease consider- 
ably. In other words, the decrease in the period is more sensitive to A than the 
temperature itself. The effect of thermal conduction on the periods compared to 
the stratified isothermal loops (case 1, A — > 1) is about 18 min. For A = 2 the 
increase in the periods is about 8 min. The deviation in the period ratio from 
unity increases with increasing L and A (for 2b a se = 6 MK). The deviation of 
the damping quality from unity in isothermal loops decreases with increasing 
L for short loops (Figure 7). However, if the loop length is further increased, 
the damping quality of non-isothermal loops becomes larger than that of the 
isothermal case. 

There is a negative relationship between the damping time and temperature in 
the case of non-isothermal loops: An increase in the temperature brings about a 
decrease in the damping time. On the other hand, there is a positive relationship 
between the damping time and the loop length. Wc sec that, in the presence of 
thermal conduction, the damping of both isothermal and non-isothermal long 
loops are in the strong regime {r&/p < 1). 

3.5. Case 5: d ^ 0, e ^ 0, r ^ 0, JV^O 

In the most general case which includes heat conduction, compressive viscosity, 
and optically thin radiation, Equation (8) must be solved. The numerical results 
are presented in Figures 8 and 9. In Figure 8, periods and their ratios are plotted 
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versus L. Compared with the case of the adiabatic loops (Figure 1), the periods 
increase noticeably. Periods with Tb aS c = 6 and 10 MK are 2-50 min and 1-45 
min, respectively (for isothermal loops). By increasing A from 1 to 2, the periods 
decrease about 20 min. 

In the presence of all dissipation mechanisms and temperature gradient, the 
damping time and damping quality are consistent with the observed values 
(Banerjee et ai, 2007; Kliem et ai, 2002; Wang et al, 2002a, 2002b, 2003, 2005). 

4. Summary 

The coronal loop has been modeled as a half-circle magnetic flux tube with a 
uniform magnetic field along the loop axis. The effect of the inhomogeneity 
of the temperature on the periods, the period ratio, and damping times of 
standing slow modes in the gravitationally stratified coronal loops was investi- 
gated in this study. The effects of optically thin radiation, compressive viscosity, 
and thermal conduction were considered. The linearized one-dimensional MHD 
equations were reduced to a fourth-order ordinary differential equation for the 
perturbed velocity. Under the solar coronal conditions (low-/3 plasma) the resul- 
tant equation was solved numerically. Moreover, the oscillations and the damping 
of various flux tube models in the presence of dissipation mechanisms were 
compared with the oscillations and damping of adiabatic and isothermal loops. 
The following shows the main results: 

• In a gravitationally stratified loop, the periods and their ratio (fundamental 
vs. first overtone) are sensitive functions of the temperature inhomogeneity. 

• In hot corona loops (T > 6 MK), the effect of optically thin radiation on 
the periods and damping of the oscillations is negligible for both isothermal 
and non-isothermal cases. 

• In the presence of compressive viscosity, the damping times significantly 
change but viscous dissipation is more effective in shorter and hotter loops. 
Compressive viscosity introduces a cutoff frequency on the oscillations. 

In the case of non-isothermal and hotter loops the computed damping times 
agree with the observed data. 

• Thermal conduction changes the periods more significantly than other dis- 
sipation mechanisms, and the damping of oscillations is in the weak regime 
(Td/p > 2). We noted that thermal conduction is a function of the tem- 
perature gradient and depends on higher derivatives, and this is the reason 
for bringing more changes in the periods in addition to the changes due to 
inhomogeneity in the temperature. 

• In the most general case where optically thin radiation, compressive viscos- 
ity, and thermal conduction are all considered, the periods are positively 
correlated with the loop length and negatively correlated with the temper- 
ature. The damping time and damping quality are complex functions of 
both loop length L and inhomogeneity parameter A. We also found that 
in the isothermal gravitationally stratified loops, thermal conduction and 
compressive viscosity are needed to reproduce the observed periods and 
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damping times. In the presence of the temperature gradient in the gravita- 
tionally stratified non-isothermal loops, the periods and the damping times 
resulted by compressive viscosity and/or thermal conduction dissipation, 
are consistent with the observed data in special cases. However, by tuning 
the dissipation parameters, the periods and the damping times could be 
made in good agreement with the observations in more general cases. 

Acknowledgements The authors thank the anonymous referee for very helpful comments 
and suggestions. 



Appendix 

Coefficient A, B, C, T>, and £ in Equation (8) are: 

T (z) = [1 + (A — 1) cos(f ], p (z) = T (z)-^ +1 \ p (z) = T Q {zf^ , 

g = sin(^), fc || = T 5/2 (z), 770(2:) = T Q 5/2 (z), 

fM--v* „[ TaP 'o I 2T «Po 7T aPo 1 35T o 2 1 7T o 1 1 r~n-T~ 2 

A(r,e,d,z)= 7 dk oll [^+T ] 

B(r, e, d, z) = 7*o|| [(^ - 2^ - + ^u] 

+ jdk oll [-Ng +^ + 2 -¥-(T + ^u)j£^] 

Civ r d "\ — I -:dh ,,rr 35T ° 2 I 7T « 14T nPo 2T "?o I 4T °p'° W f 4 7-1 

+ 7&o||(5-^-|)W4^-^o] 

+ 7*o||[^ + ^ + ^ + ^-^-2iV 5 '-«-7, 2 ] 
+ 7 rffeo||[iV ff - 77^ - ^ + (-^ + )7%4- ~ ^^jflfa 
+ 3r 7/?u T - 1 + (/(z) - i/)^!/ - 



(r, e, d, z) = +7*o|| [§£ + ^ " T " T + T X^f" - ^gpo)] 

1U ..\ 21T oPo i ST 0Po' 9T »PoPn 21T nPo 2 i 6TqPq 3 . 21p' Tg 

■^ dk o\\[—^ + — ^ p-^ + ^T + ^- 



4T 4p T T 2 t gT( 

' U N L 3po » p Po Po 1 

+ 7*o|| [(^gf - ^ - ^)(7^^ - A^Vo - 2iV. 9 p - 1Po v 2 )] 
+ Sr^KSpoPoTo- 1 - ^T^T - 2 ) - 

- (1 + 7)Po" + 7Po^ rjgy^ + (/(*) - Ngp ) 
£(r,e,d,z) = -7*0111(5 " ^ " ^) W + Wo + 7^ 2 ) 

- 7*o|| [S + ? - ^ - ^ + ^][« + 7Po, 2 ] 

-L^Jh r 63T oPoPo I 21T oPo 3 i mST^p^ 12T Q pgp'£ 6T p^ 4 1057^ lOBT^pg 
+ 7*>||l 27 5 r _ + ^^+ 8p 7- + p| Ip^T+ 4p T„ 
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_ 3TqPq 2 _ iTop'aPg' 2lTg P {" Pq"T _ 21 p' a Tg p' 105T Tp' p' 7Tg" p' B 21T n , p ' i 
p'f ~~pj r 2p ^ po 2p^ r 4T po x 2p ^ 2p J 

-7*o||(^ + ^ + ^ + ^ 2 ) 

+ 7*0,1 + ^ + > 2 + ~ 2 ^)(N 9 p' + 7A^)] T7 fc 

+ Sry^ 2 ^ 1 - p Po^T - 2 + poTq-Vo') - %%?- f X Y ) 
- ^o+Po^xm^J ~ " p K N 9Po + 7Po^ 2 ) 

Here all quantities are dimensionless. The prime indicates derivatives with re- 
spect to z. 
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Figure 1. The periods of the fundamental mode pi (top row) and of the first overtone p2 
(middle row), and their ratios (bottom row) for case 1 (see text) are plotted versus L in 
oscillating stratified loops with different values of A and ?b ase . 
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Figure 2. The periods of the fundamental mode pi (top row) and of the first overtone 
P2 (middle row), and their ratios (bottom row) for case 2 (see text) are plotted versus L in 
oscillating loops in the presence of radiation and gravity for different value of A and 7b ase . 
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Figure 3. The damping times of the fundamental mode (rdi) and of the first overtone (t<j2) 
(top row), and the damping times per period {r^i/pi, i"d2/P2) (bottom row) for case 2 (see 
text) arc plotted versus L in oscillating loops in the presence of radiation and gravity for 
different value of A and 7b asc . 
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Figure 4. The periods of the fundamental mode p\ (top row) and of the first overtone 
P2 (middle row), and their ratios (bottom row) for case 3 (see text) are plotted versus L in 
oscillating loops in the presence of viscosity and gravity for different value of A and 7t, ase . 
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Figure 5. The damping times of the fundamental mode (rdi) and of the first overtone (t^) 
(top row), and the damping times per period (t^i/pi, T < j2/p2) (bottom row) for case 3 (see 
text) are plotted versus L in oscillating loop in the presence of viscosity and gravity for different 
value of A and Tb asc . 
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Figure 6. The periods of the fundamental mode pi (top row) and of the first overtone 
P2 (middle row), and their ratios (bottom row) for case 4 (see text) are plotted versus L in 
oscillating loops in the presence of conduction and gravity for different value of A and 7b asc . 
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Figure 7. The damping times of the fundamental mode (rdi) and of the first overtone (t^) 
(top row), and the damping times per period (fdi/pi , id2/:P2) (bottom row) for case 4 (see 
text) are plotted versus L in oscillating loops in the presence of thermal conduction and gravity 
for different value of A and 7b ase . 
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Figure 8. The periods of the fundamental mode pi (top row) and of the first overtone 
P2 (middle row), and their ratios (bottom row) for case 5 (see text) are plotted versus L in 
oscillating loops in the presence of radiation, viscosity, thermal conduction, and gravity for 
different value of A and 7b asc . 
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Figure 9. The damping times of the fundamental mode (t<u) and of the first overtone (t^) 
(top row), and the damping times per period (r^i/pi, t^/vz) (bottom row) for case 5 (see 
text) are plotted versus L in oscillating loops in the presence of radiation, viscosity, thermal 
conduction, and gravity for different value of A and T basc . 
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